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Bianchi type I magnetized cosmological models in the presence of a bulk viscous fluid are investigated. 
The source of the magnetic field is due to an electric current produced along x-axis. The distribution 
consists of an electrically neutral viscous fluid with an infinite electrical conductivity. The coefficient 
of bulk viscosity is assumed to be a power function of mass density. The cosmological constant A is 
found to be positive and is a decreasing function of time which is supported by results from recent 
supemovae observations. The behaviour of the models in presence and absence of magnetic field are 
also discussed. 

1. Introduction 

The occurrence of magnetic fields on galactic scale is well-established fact today, and 
their importance for a variety of astrophysical phenomena is generally acknowledged as 
pointed out Zeldovich et al.13 Also Harrison □ has suggested that magnetic field could have 
a cosmological origin. As a natural consequences we should include magnetic fields in the 
energy-momentum tensor of the early universe. The choice of anisotropic cosmological 
models in Einstein system of field equations leads to the cosmological models more general 
than Robertson- Walker model.Q The presence of primordial magnetic fields in the early 
stages of the evolution of the universe has been discussed by several authors. ~0 Strong 
magnetic fields can be created due to adiabatic compression in clusters of galaxies. Large- 
scale magnetic fields give rise to anisotropics in the universe. The anisotropic pressure 
created by the magnetic fields dominates the evolution of the shear anisotropy and it decays 
slower than if the pressure was isotropic. QEl Such fields can be generated at the end of 
an inflationary epoch. Anisotropic magnetic field models have significant contribu- 

tion in the evolution of galaxies and stellar objects. Bali and AliEll obtained a magnetized 
cylindrically symmetricjinhierse with an electrically neutral perfect fluid as the source of 
matter. Several authors t3~ta have investigated Bianchi type I cosmological models with a 
magnetic field in different context. 
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Most cosmoiogicai models assume that the matter in the universe can be described by 
'dust' (a pressure-less distribution) or at best a perfect fluid. Nevertheless, there is good 
reasonto believe that - at least at the early stages of the universe - viscous effects do play 
a role.E2l^cJ For example, the existence of the bulk viscosity is equivalent to slow process 
of restoring equilibrium states.El The observed physical phenomena such as the large en- 
tropy per baryon and remarkable degree of isotropy of the cosmic microwave background 
radiation suggest analysis of dissipative effects in cosmology. Bulk viscosity is associated 
with the GUT phase transition and string creation. Thus, we should consider the presence 
of a material distribution other than a perfect fluid to have realistic cosmoiogicai models 
(see Gr0r£j for a review on cosmoiogicai models with bulk viscosity). The effect of bulk 
viscosity on the cosmoiogicai evolution has been investigated by a number of authors in the 
framework of general theory of relativity 

Models with a relic cosmoiogicai constant A have received considerable attention re- 
cently among researchers for various reasons (see Refs. 0-0 and references therein). Some 
of the recent discussions on the cosmoiogicai constant "problem" and on cosmology with 
a time-varying cosmoiogicai constant by Ratra and Peebles£^ Dolgov@-il and Sahni and 
StarobinskyU point out that in the absence of any interaction with matter or radiation, the 
cosmoiogicai constant remains a "constant", however, in the presence of interactions with 
matter or radiation, a solution of Einstein equations and the assumed equation of covariant 
conservation of stress-energy with a time-varying A can be found. For these solutions, 
conservation of energy requires decrease in the energy density of the vacuum component 
to be compensated by a corresponding increase in the energvjlensity of maltei or radiation. 
Earlier researchers on this topic, ^e contained in ZeldovichJSj Bertolami, La€ZlWeinberg,EJ 
and Carroll, Press and Turner. Recent observations by Perlmutter et al. i3 and Riess 
et al. Ej strongly favour a significant and positive A. Their finding arise from the study 
of more than 50 type la supernovae with redshifts in the range 0.10 < z < 0.83 and 
suggest Friedmann models with negative pressure matter such rs a cosmoiogicai constant, 
domain waUs or cosmic strings (Vilenkin, El3 Garnavich et a/.Ell). Recently, Carmeli and 
KuzmenkoEd have shown that the cosmoiogicai relativity theory (Behar and CarmeliEJ) 
predicts the value A = 1 .934 x lO^^^s^^ for the cosmoiogicai constant. This value of A is 
in excellent agreement with the measurements recently obtained by the High-Z Supernova 
Team and Supernova Cosmoiogicai Project (Garnavich et a/.;Ell Perlmutter et a/.;Ej Riess 
et flZ.;E3 Schmidt et alM). The main conclusion of these works is that the expansion of the 
universe is accelerating. 

Several ansatz have been proposed in which the A term decays with time (see Refs. 
GasDerini£jii3 Berman,EJ Freese et ai. @ Ozer and Taha@ Peebles and Ratra,E3Chen and 
Hu,E3 Abdussattar and Viswakarma,Ll] Gariel and Le Denmat,Lil Pradhan et al.\ m^). Of 
the special interest is the ansatz A cx S^^ (where S is the scale factor of the Robertson- 
Walker metric) by Chen and Wu,EI which has been considered/modified by several authors 
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( Abdel-Rahaman^ Carvalho etal., @ Wagai Silveira and Wagafi Pradhan and Pandeyfl 
VishwakarmaLj) . 

Recently Bali and GokhrooEl obtained a Bianchi type I anisotropic magnetized cosmo- 
logical model for perfect fluid distribution. Motivated by the situations discussed above, in 
this paper, we shall focus upon the problem of establishing a formalism for studying the 
general relativistic evolution of magnetic homogeneities in presence of bulk viscous in an 
expanding universe. We do this by extending the work of Bali and Gokhroo0 by including 
an electrically neutral bulk viscous fluid as the source of matter in the energy-momentum 
tensor. This paper is organised as follows. The metric and the field equations are presented 
in Sec. 2. The Sec. 3 includes the solution of the field equations and physical and geometrical 
features of the models. In Sec. 4, we obtain the bulk viscous cosmological solutions in the 
absence of the magnetic field. In Sec. 5, we discuss our main results and summarize the 
conclusions. 



2. Field Equations 

We consider the Bianchi type I metric in the form 

ds^ ^ -dt^ + A^dx^ + B^dy^ + C^dz^, 



(1) 



where A, B and C are functions of t only. The energy momentum tensor in the presence of 
bulk stress has the form 

Tl = {p+p)v,v^ +Pgl+El (2) 
where is the electro-magnetic field given by LichnerowiczH as 



and 



H ( v.v^ + -9l 



- h,h^ 



(3) 



(4) 



Here p, p, p and ^ are the energy density, isotropic pressure , effective pressure, bulk viscous 
coefficient respectively and is the flow vector satisfying the relation 



gijv'v^ 



-1 



/i is the magnetic permeability and hj the magnetic flux vector defined by 



/i, = i *Fj,v^ 



where *Fij is the dual electro-magnetic field tensor defined by Synge Lil to be 



0, 



9 pki 



(5) 



(6) 



(7) 
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Fij is the electro-magnetic field tensor and etjki is the Levi-Civita tensor density. Here, the 
comoving coordinates are taken to be = Q = = and = 1. We take the incident 
magnetic field to be in the direction of x-axis so that hi 7^ 0, /i2 = = = /i4. Due to 
assumption of infinite conductivity of the fluid, we get F14 = = F24 — F34. The only 
non- vanishing component of Fij is F23. The first set of Maxwell's equation 

Fij-k + Fjk-i + Fki-j = (8) 

leads to 

F23 = /(const.) (9) 
where the semicolon represents a covariant differentiation. Hence 



AT 

- JBC '''' 



The Einstein's field equations 



Ri - ^Rgi + ^gi = -SttT/ , (c = l, G = l in gravitational unit) (11) 
for the line element ( 1 ) has been set up as 



in[p~ ^_ ) ^ ^-^-^^-A (12) 

(13) 



/2 ^ 


, B44 


C44 


B4C4 




'^B 


^ C ^ 


BC 


f ^ 


, A44 


C44 


A4C4 








AC 


- \ 


A44 


B44 


A4B4 


IjlB^C^ ) 




" ~B~ " 


AB 




A4B4 


A4C4 
\ - 


B4C4 
1 ^— ^ 


IfiB^C^ J 


AB 


AC 


BC 



^[P+ = r-^-^TT-^ (14) 



The suffix 4 by the symbols A, B and C denote differentiation with respect to t. 
3. Solution of the field equations 

Equations ([l^) - ( p^ represent a system of four equations in six unknowns A, B, C, p, 
p and A. To get a determinate solution, we need two extra conditions. First we assume 

A^BC (16) 

From Eqs. ([ij) and ([l4|), we get 



B44 C44 A4 I C4 B4 

~B c" ^ T \~c ^ "b 

which leads to 

1/4 a 
V e 



(18) 
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where BC ^ e, ^ = v, and a is a constant of integration. From Eqs. ( p^ ) and (pj[), we get 

K ^44 -^44 A.4C4 I34C4 



B^C^ A B AC BC 

where 



By using the condition ( jl^ in ( 19), we have 

C44 IB4C4 C? iiT 



C BC C2 52(72 

which leads to 

£44 el _K 
2e 2e2 " 
Setting £4 = /(e) in Eq. (0), we have 



Integrating Eq. ( p4| ) leads to 



From Eq. dlq), we get 

V4 a 



V {afi + bt + d) 

Integrating Eq. (^) leads to 
1 



where M is a constant of integration. 



(19) 



i. = -M! (20) 



(21) 



(22) 



2eff + 2f = 4K (23) 

which reduces to 



+ (24) 



where N is the constant of integration. Eq. ( |2q ) on integration leads to 

v/(2-ft:e2 + iV) = 2Kt + L, (26) 
where L is the integrating constant. Eq. (EN) can be rewritten as 



e = V(ai^ + 6i + rf) , (27) 

where ^ 

a = 2ii:; 6 = 2L and d = i^L_iXl. (28) 

2K ^ ' 



(29) 



= — eV(4»d-i.^) V(4a<i-b2) if4ad-b^>Q (30) 
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where ^ > 0. 



if 4ad -b^ = 



3.1. Case (i) : when 4ad — b^>0. 



In this case, we have 



^at^ + bt + d 



,-1 (2at + b) 



J (at^ + bt + d) , ^° , tan- . 



M 



M^{afi + bt + d) 



■ 2c, (2at + b) 



Hence, the geometry of the spacetime (|l|) reduces to the form 

ds^ ^ -dt^ + [at^ + bt + d)dx^ 
yj{at^ + bt + d) 



M 



M y/ {at^ + bt + d) 



• ^(iad-b-) 



-I (2at + 6) 
\/(4ad-62) 



The effective pressure p and density p for the model ( pq ) are given by 



87rp = 



16{at^ + bt + d)^ 
5(2at + b)^ -4a^ 



2{at^ + bt + d) 



-A 



K 



I6{afi + bt + dy 2{afi + bt + d) 
where 9 is the scalar of expansion calculated for the flow vector i;* and is given by 

g _ {2at + 6) 
{afi + bt + d) 



(32) 



(33) 
(34) 

(35) 



(36) 



(37) 
(38) 

(39) 



For the specification of ^, now we assume that the fluid obeys an equation of state of the 
form 

P = IP, (40) 

where 7(0 < 7 < 1) is a constant. 

Thus, given £^{t) we can solve the cosmological parameters. In most of the investigarions 
involving bulk viscosity is assumed to be a simple power function of the energy density.E3~a 



m = Cop", 



(41) 
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where and n are constants. I£n = 1, Eq. (^ij) may correspond to a radiative fluidH 
However, more realistic modelsLj are based on lying in the regime < n < 3. On usmg 
Eq. @ in Eq. @, we obtain 



where T ^ af + bt + d. 
3.1.1. Model (i) : = Co)- 



(42) 



When n = 0, Eq. reduces to ^ = ^o- With the use of Eqs. (|38|), (|39|) and Eq. 
( |4^ reduces to 



87r7p = 



^7r(2at + 6)^0 4(&2 - 4ad) + (2at + bf 



T 



16T2 



2r 



-A 



Eliminating p{t) between Eqs. ( |3^ and (|43|), we get 

1 



(1+7)A: 



16T2 

1 r 

3.1.2. Mode? fii; ; = ^op) 



(2at + 6)2(1 - 57) + 4(6^ - 4afi) + 470^ 
K{1 -7) + 16<()(2ai + &) 



(43) 



(44) 



When n = 1, Eq. (|41|) reduces to ^ = ^op. Eq. (^, with the use of Eqs. (|38|), (g9|) 
and (p^, reduces to 



{2at + b)^o 



rp2 



87rp 



Eliminating p{t) between Eqs. ( |38| ) and (|45|), we get 

1 



4(&2 - 4ad) + (2a< + bf - 4a2 



16^2 



2r 



(45) 



16 [(1 + 7)^2 -(2ai + 6)^0] 
Aa^ - 8KT - {jT^ - (2at + b)^o 

\ 5 (2at + bf - 4a^ + 8KT 
X 



4 (5^ - 4-ad) + (lat + bf 



(46) 



From Eqs. (44) and (|4^, we observe that the cosmological constant is a decreasing function 
of time and it approaches a small value as time progresses (i. e., the present epoch), which 
explains the small value of A at present. In the model (i), we also find that A is positive for 
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all values of 7 lying in the interval (0, ^) and > 5{4ad — b^). In the model (ii), we see 
that the second term of the Eq. (^6|) decreases rapidly with expansion of the universe and 
the first term will dominate over the second term. Thus, from Eq. (^6|) it is also observed 
that A is positive when K < ^ and > ^Kd + 4a^. These small positive values of A 
in both models (i) and (ii) are supported by results from recent supernovae observations 
(Perlmutter et al., [58] Riess et al., [59] Garnavich et al., [61] Schmidt et al. [64]). 
Physical and Geometrical Features of the Models 
The component of the shear tensor cr^ are given by 

I _ {2at + b) 



6T 




{2at + b) 


a 


12T ^ 


2T 


{2at + b) 


a 


12T 


2T 


at = 





Thus 



(47) 



, {2at + bf ~ 12^2 

cr^ = -^^ ^-^ (48) 

48r ^ ' 



The rate of expansion Hi in the direction of x, y, z -axes are given by 

^ ^ (2at + &) 



2T 




{2at + b) 


a 


AT ^ 


^ 2T 


{2at + b) 


a 


AT 


" 2T 



H3 = ^ — (49) 

4r 2T 

The model starts expanding att > — ^ and the expansion in the model decreases as 
time increases and stops at t = 00. The model, in general, represents shearing and non- 
rotating universe. Since limt^oo f 7^ 0, therefore, the model does not approach isotropy 
for large values of t. There is a singularity in the model at t = which is real physical 
singularity. 

3.2. Case (ii) : when Aad- b^ < i.e. Aad- b^ ^ -I, I > 0. 



In this case the geometry of the spacetime (|I]) takes the form 
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1 

2 


"(^^ 


- -^1 

2a) 


_ VI' 

2a 






2a) 


^ 2a J 



The effective pressure and density for the model ( pO| ) are given by 

4a2t2 + 5&2 - 4a6i - I6ad K 



87rp: 



5{t 



_b_\2 

2a) 



16Ti 2^aTi 



-A 



4T2 



16Ti 2/iaTi 



where the scalar expansion is given by 



and 



T, = 



2ft + ^) 



(50) 

(51) 
(52) 

(53) 



3.2.1. Model (i): (C = ^o). 



When n = 0, Eq. (^ij) reduces to ^ = and hence Eq. (|TJ), with the use of Eqs. (p^, 
|) and leades to 



AaH^{l + 57) + 56^(1 + 7) - 4at(l - 567) - 16arf 



167r7p = 167r^o(^ + tt- 



2a \6a?T'^ 



^1-7) if(l+7) 



\6T\ IfiaTi 
Eliminating p{t) between Eqs. ( |5^ ) and (^J), we get 

1 



- (1 - 7)A 



(54) 



(1+7)A 



4a¥(l - 57) + 55^(1 - 7) - 4a6t(l + 57) - I6ad 
a2(l-7) , if(l-7) 167reoft+^) 



3.2.2. Model (ii) 



\6T\ 2/xaTi 



T, 



(55) 



When n = 1, Eq. (|lj) reduces to ^ = ^oP- The Eq. (|T]), with the use of the Eqs. @, 
3) and (p3|), reduces to 



87rp[a(l +7)ri - (2at + 6)^0] 



Aa^f + 56^ - 4a6t - 16ad 



16ari 



5(2ai + 6)2 aa^ K 



(56) 
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Eliminating p{t) between Eqs 

A[a(l+7)Ti-(2ai + fe)eo] 



and (p4|), we get 

AaH'^ + 56^ - Aaht - \6ad aa^ K 



16aT, 



16 2/z 



a'yTi — {2at + 



4T2 



2/iaTi 



(57) 



In both models (i) and (ii), we have observed from Eqs. ( p5[ ) and ( p7[ ) that the cosmological 

constant is a decreasing function of time and it approaches a small value as time progresses 

(i.e., the present epoch), which explains the small value of A at present. 

The Physical and Geometrical Features of the Models 

The non-vanishing components of shear tensor (cr^ ) are given by 



1 

12 
1 

"l2 



3a 2(t+A) 



I 

4a? 



3a 



+ 



(58) 



(59) 



(60) 



The models start expanding with a big bang at i = — ^ + 2^ and the expansion in the models 
decreases with time. The expansion in the models stops at t = 00. Since limt^oo^ 7^ 0. 
Therefore, the models do not approach isotropy for large values of t. 
3.3. Case (Hi) : when 4ad ~ — . 



In this case the metric (fil) reduces to 



where 



ds'- = -dT^ + T| dX^ + e dY^ + T2 



T2^t+ — 

2a 



The effective pressure p and density p for the model (^Tj) are given by 



K 



4T| Aa^T^ 2aTl 



K 



4r| Aa-T^ 2aTl 



where the scalar expansion 9 is given by 



2 



(61) 



(62) 
(63) 

(64) 
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3.3.1. Model (i): = 6 



When n = 0, Eq. ( |41| ) gives ^ = and hence Eq. 
and (|64|), reduces to 



1 



^ 877(1 +7)T2 
Eliminating p{t) between Eqs. ( |63| ) and (|6^, we get 

16^0 



2T2 2a^Ti 



with the use of Eqs. (M,, (b3b 



(65) 



(1+7)A: 



To 



12 ~ru,.2 

3.3.2. Mode/ (ii) : = 



(66) 



When n = 1, Eq. (0) gives ^ = ^oP and hence Eq. (||), with the use of Eqs. (0), (l^]) 
and (^J), reduces to 



a' 



2 1 



8^[(l+7)T2-2^o] 
Eliminating p(t) between Eqs. ( |63| ) and (|67|), we get 



2T2 2a2T,3 



(67) 



A = 



{(l+7)72-2Co}T2 



(1-57) (1+7)^ (l-7)a^ 



2a 



4a2T2 



2^0 r 5 X 



r2 [4 2a 



(68) 



From Eqs. ( |66[ ) and (|68|), we observe that the cosmological constant is a decreasing function 
of time and it approaches a small value as time progresses (i. e., the present epoch), which 
explains the small value of A at present. In the model (i), from Eq. (|6|), it can be seen 
that first and second terms dominate over the third term. Further we also find that A is 



a(l-57) 
2(1+7) 



In the model (ii). 



positive for all values of 7 lying in the interval (0, i) and K < 

from Eq. (|6^), one can easily see that that for all value of 7 lying in the interval (0, i) and 

^0 < ■^^^^^ K < °'2l-i+^) ' A positive. These small positive values of A are consistent 

by the results from recent supernovae observations. 

The Physical and Geometrical Features of the Models 

The non-vanishing components of shear tensor (cr^ ) are given by 



1 



3a 



1 

3a 1 



(69) 
(70) 
(71) 
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The models start expanding with a big bang t = and the expansion in the models 
decreases with time increase. The expansion in the models stop at i = cx3. Since Hmt _,oof- 7^ 
0. Therefore, the models do not approach isotropy for large values of t. 
4. Solution of the field equations in absence of magnetic field 



In the absence of magnetic field, the metric (^6|) reduces to 



ds^ = -dt^ + ibt + P)d.^ + W^dy^ + ^[VE±Rd^\ ill) 

where 

Thus when ii' ^ then d ^ (3. In absence of the magnetic field i.e. when K ^ Q then 
effective pressure p and density p for the model ( [72| ) are given by 

8., = 84, - = + ^^^^ - - A, (74) 



^2 

16{bt + 4{bt + P) 

where scalar expansion 6 is given by 

b 



8-^=7771^-^- 7777— A, (75) 



(76) 



bt + (3 

4.1. Model (i): (C = Co)- 



When 71 = 0, Eq. (|4ip gives ,^ = and hence Eq. ([74|), with the use of Eqs. (^, (|75| 
and (|76|), reduces to 

^7r6Co - - 2d 

{bt + 13)^ ^bt + P) 

Eliminating p{t) between Eqs. ( [t^ and ([77|), we get 



8^(1 +7)P=7ir^ + ^71^-77^ (77) 



4.2. Model (ii) : (C = Cop)- 



When n = 1, Eq. (|T]) reduces to ^ = ^op. Eq. @, with the use of Eqs. @, @, 
and (|76|), reduces to 

- (5^2 ^ 4q2 - 2d) 

^""P- 8{bt + p)[{\+^){bt + P)-b^o] ^ ' 
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Eliminating p{t) between Eqs. ( |7^ and (|7^, we get 

(1 -7)(562-4a2) _4a(i- 



(5b--4a-)&^o 
ibtW) 



\6{bt + (3) [(l+7)(6t + /3)-6eo] 



(80) 



In both the models (i) and (ii), from Eqs. ( |78| ) and (pOj), we observe that the cosmological 
constant is a decreasing function of time and it approaches a small value as time progresses 
(i. e., the present epoch), which explains the small value of A at present. In the model (i), 
from Eq. ([ts]), it can be seen that for all time t > ^^^^^^^ — ^, the cosmological "constant" 
is positive definite. In model (ii), Eq. ( |8C| ) suggest that the positivity of A demands the 
relation among constant as 6^ > | {a^ + ad) and < '-^^^''^ . These small positive values 
of A in both models are supported by the results from recent supernovae observations. 
The Physical and Geometrical Features of the Models. 
In the absence of magnetic field, the shear {af ) are given by 

Here, limt^oc f 7^ 0. The model does not approach isotropy for large values of t in absence 
of magnetic field also. 
5. Conclusions 



We have obtained a new class of Bianchi type I anisotropic magnetized cosmological 
models with a bulk viscous fluid as the source of matter. Generally, the models are 
expanding, shearing and non-rotating. In all these models, we observe that they do not 
approach isotropy for large values of time t either in the presence or in the absence of 
magnetic field. 

The cosmological constant in all models given in Sec. 3 and Sec. 4 are decreasing function 
of time and they all approach a small value as time increases (i.e., the present epoch). 
The values of cosmological "constant" for these models are found to be small and positive 
which are supported by the results from recent supernovae observations recently obtained 
by Ihe High - Z Supernova Team andSupernova Cosmological Project ( Garnavich et 
flZ.;Lll Perlmutter et fl/.;E3 Riess et aZ.;E3 Schmidt et al. Ej). Thus, with our approach, 
we obtain a physically relevant decay law for the cosmological constant unlike other 
investigators where adhoc laws were used to arrive at a mathematical expressions for the 
decaying vacuum energy. Thus our models are more general than those studied earlier 
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